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Abstract 

The construction of ready to use templates for gravitational waves from spinning binaries is an important 
challenge in the investigation of detectable gravitational wave signals. Here we present a method to evaluate 
the gravitational wave polarization states for inspiralling compact binaries in the extreme mass ratio limit. 
We discuss the effects caused by the rotation of the central massive object for eccentric orbits in the Lense- 
Thirring approximation and give the formal expressions of the polarization states including higher order 
corrections. Our results are in agreement with existing calculations for the spinless and circular orbit limits. 

1 Introduction 

Compact objects forming binary systems are among the possible sources of gravitational radiation, whose 
detection is expected by gravitational wave observatories, i.e. LIGO [T], VIRGO 2 , GEO600 3], TAMA [3] 
and LISA |5| . Under radiation reaction the orbit of the binary slowly decays and the frequency and amplitude 
of the emitted signal increase over time. For the evaluation of the the noisy output of the detector, the method 
of matched filtering [6] is employed, whose effectiveness depends on the precise knowledge of the gravitational 
waveform. 

One of the commonly used approximation schemes for the computation of the dynamics and the emitted 
gravitational wave signal is the post-Newtonian (PN) expansion [71 [51 IS]- In this weak field approximation, the 
velocities and the gravitational potential involved are small, but there is no restriction on the mass ratio of the 
components. This approach is considered to be a precise tool to describe the motion of an inspiralling binary up 
to the innermost stable circular orbit, which is determined with high accuracy in [10]. Close to the merger phase 
the PN approximation becomes less accurate, known as the intermediate binary black hole problem [111 I12j. 
therefore one has to use different techniques. 

The computation of the gravitational wave polarization states was carried out e.g. in [131 1141 I15j and 
[16[ I17[ riSj for quasi-circular and elliptic orbits, respectively. For spinning binaries the evaluation of the wave 
pattern has been discussed by several authors, see e.g. [9l[19l[20j. In most cases the gravitational waveform is 
expressed formally, i.e. in terms of the dynamical quantities of the motion. Up to now different calculations 
which have determined the explicit time dependence of the emitted waves analytically are restricted to the 
circular orbit or spinless cases. 

Here we present a method to express the gravitational waveforms of a binary system using the PN approx- 
imation. The motion of the binary is described in the Lense-Thirring approach [21^, when we consider the 
geodesic motion around a central spinning body. We give the formal expressions of the wave polarization states 

and up to 1.5 relative PN order beyond the leading order Newtonian expressions in the case of eccentric 
orbits. To describe the effects of rotation we focus on terms linear in the spin of the central, massive body. 
Moreover, we reproduce the circular orbit case in order to determine the influence of the eccentricity on the 
detectable gravitational waves. 

Our work is based on the results of Kidder, Will and Wiseman [HI US]- With the use of the general definition 
of the mass and current multipole moments they have determined the form of the transverse-traceless tensor 
hl^rp representing metric perturbations up to 1.5 relative PN order in the case of a spinning compact binary in 
terms of the constants of the motion, the separation and relative velocity vectors and the vector representing the 
direction of the line of sight. Since the expressions contain these quantities in a formal way, their results can be 
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used both in the eccentric and circular orbit cases. They have given a very useful form of the transverse-traceless 
tensor in the circular orbit case with the decomposition of the relative velocity vector into components parallel 
and perpendicular to the separation vector. We give a similar decomposition of hlfj, in the eccentric orbit case 
and we use their results to check our formulae. 

To discuss the motion of the orbiting bodies we use the results of Gergely et. al. [22], where the complete 
radial and angular description of the motion is given in the Lense-Thirring approximation. The Lagrangian 
formalism was chosen as the starting point in determining the equations of motion. They have given an 
appropriate parameterization of the orbit (for further details see |23j). which can be the basis of the integration 
of the dynamics. 

In order to avoid long and complicated expressions we use the comoving coordinate system, where the x axis 
is aligned with the separation vector r and the z axis with the newtonian angular momentum Ln = /xr x v. Here 

V is the relative velocity vector. We describe the motion in the comoving system since the transverse-traceless 
tensor, /i^^,, and the projection formulae are simple in this case. For a convenient description of the motion 
we introduce a coordinate system which does not change over time, called invariant system. We can fix the z 
axis of this coordinate system to the spin vector of the massive body, since the precession of the spin can be 
neglected in the Lense-Thirring approximation. 

We organize our paper as follows. In Sec. II we describe the dynamical elements of the motion in the 
Lense-Thirring approximation, in the invariant and comoving coordinate systems. For the description of the 
time evolution of these elements we introduce a third coordinate system where the x and z axes are aligned 
with the separation and angular momentum vectors, r and L. The transformation between the invariant and 
comoving systems is parametrized with Euler-angles. In Sec. Ill we compute the explicit form of the vectors 
which describe the relative position of the binary and the detector in order to express the polarization states. 
Moreover we give the components of the spin vector in the comoving system. Sec. IV contains our main results, 
namely the splitting of the detectable gravitational wave signal into two polarization states /i+ and and 
their derivation from /i^^n. To complete the general description of the method we show explicit contributions 
to the gravitational wave signal which belong to the different polarizations, PN orders and spin effects. In Sees. 

V and VI we discuss the spinless and circular orbit limits of our calculations. 

Throughout the paper we use the c = G = 1 convention. 



2 Description of the motion 

The components of the hl^j, tensor have the simplest form in a comoving system, where the x axis lies in the 
direction of the separation vector r and the z axis in the direction of the Newtonian angular momentum Ln. 
Since the polarization states are scalar quantities, we can choose this comoving system to determine their form. 

In the Lense-Thirring approximation let M and S denote the mass and the spin vector of the rotating central 
body and m the mass of the particle (m ^ M) orbiting it. 

When the two bodies have comparable masses the evolution of the spin is governed by the spin precession 
equations [25] . Since the length of the spin vector S is constant the relative order of the change of the spin is 

S--e, (1) 

where e is the PN parameter. When the mass ratio is negligible the precession of the spin vanishes. In this case 
the z and y axes of the invariant system can be aligned with the direction of S and S x N, respectively, where 
N denotes the direction of the line of sight. 

To describe the transformation between the comoving and the invariant systems we use Euler angles [221 124j . 
The vector u in the comoving system can be parametrized in the invariant system as 

u' = i?,($)i?,(/,jv)i?.Wu , (2) 

where (-at is the angle between S and Ln, represents the precession of the orbit (precession of Lat over S), 
and \1/ is the polar angle in the orbital plane. In the invariant system the separation vector has the form 

cos $ cos 5* — cos sin <I> sin '5 

sin $ cos 5* -I- cos l n cos $ sin ^E* | . (3) 
sin sin 5* 

Using the notation of [22j the radial equation of motion can be written as 

2 2E 2M AL,S 
r = — + FT — ' 4 
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where the total energy the magnitude L = |L| and the z component Lz of the angular momentum L = 
Ljv + L50 are constants of the motion (the expression for the spin-orbit part L50 of L can be found in [22] ). 

In order to get the equations for the angles, we introduce a third coordinate system called secondary system, 
similar to the comoving one. In this system the x axis is aligned with the separation vector r and z with the 
angular momentum L. Again, every vector u of this new coordinate system in the invariant one become 



(5) 



Here l' is the angle between S and L. Since the magnitude and the z component of L {L and respectively) 
are constants of the motion, l' is constant up to 1.5 PN order. The angle $' corresponds to the precession of L 
over S, and 5*' shows the direction of r in the plane perpendicular to L. The equations of the motion for these 
angles are given as 



L 



2S 



(6) 



The solutions of Eq. 
3): 



determine the equations of motion for the original angles (for further details see 



4- = 



L 



S cos i' 
Lr3 



{2L - mrf sin2\l>') 



25* sin*' 
25* sin i' cos* 



Lr3 



[mrf cos *' + 2L sin 

(mrf cos *' + 2L sin . 



For later convenience we decompose the relative velocity vector as 

V = f n + tij_m , 



(7) 



(8) 



where n = r/r is the direction of the x axis of the comoving system and m denotes the direction of its y axis. 
The relative velocity vector in this system is 



r(cos tAr$ 




(9) 



From the expressions of the constants of motion we can determine the magnitude and the components of the 
relative velocity vector (see [M]). Using the decomposition of the length of the separation vector r = ro 
these components are 



Wo 



Mri 



ir)o 



Mrim?rQ — L^ri + 2Lzm 



S , 



Lmrff 



S 



where we have introduced the short-hand notations 

f2E 2M 
Wo = \/ 1 , {r)o = 




riS 
(10) 



(11) 



3 Basic vectors in the comoving system 

To describe the projections of hl^j, for the calculation of the polarization states h+{t) and hx(t), we determine 
the components of the orthonormal triad (N, p, q). The vector N denotes the direction of the line of sight, 
p is chosen to be perpendicular to the Newtonian angular momentum (this way p lies in the direction of the 
node line, the intersection of the orbital plane of the source and the plane perpendicular to N called plane of 
the sky), and q = N x p. 

In the invariant system N = (sin 7, 0, cos 7), where 7 is the constant angle between N and S. This way the 
components of N in the comoving system are: 

cos * cos $ sin 7 — sin * cos l jv sin $ sin 7 + sin * sin l n cos 7 \ 

- sin * cos $ sin 7 — cos * cos iAT sin <& sin 7 -I- cos * sin tjv cos 7 . (12) 
sintjvsin$sin7 + costAf COS7 / 
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There are three conditions for the vector p. It is a unit vector which is perpendicular to N and La?. From 
the last condition we get that Pz = in the comoving system. With the use of the other conditions, the form 
of p in the comoving system becomes 



1 



sin ^E* cos $ sin 7 + cos ^ cos t n sin $ sin 7 — cos 5' sin l n cos 7 



p = — cos ^I^ cos $ sin 7 — sin cos ijv sin $ sin 7 + sin sin ijv cos 7 , (13) 

N \ g 



where 



N = JN^ + N^ ^ ^1 - = yj\ - (sinijv sin $ sin 7 + cos t at cos 7)2 . (14) 



In order to determine the polarization states, we will also need the form of the spin vector in the comoving 
system, 

(sin^'siniAT \ 
cos^sintAT . (15) 
COSiAT j 

In this way we have determined all the components of the basic vectors needed to evaluate the projections 
of the transverse-traceless tensor in terms of the elements of the motion. 



4 The formal expressions of the polarization states 

The signal h(i) of a laser-interferometric gravitational wave detector can be decomposed into two polarization 
states [20] /i+(t) and h^(€) 



h{t) - F+ia, P,Oh+{t) + F^{a, P,Oh^it) 



(16) 



The beam-pattern functions and Fx depend on the direction of the source (angles a and f3) and the 
polarization angle {£,): 



= -(l+cosa ) cos 2/3 cos 2^ + cos a sin 2/3 sin 2^ , 
Fx = — -(1 + cos a^) cos 2/3 sin 2^ + cos a sin 2/3 cos 2f . 
The independent polarization states of the waveform can be projected out from hl^rp as 



P -ft'^ 



The components of the projectors have the form: 



In the post-Newtonian approximation up to 1.5 PN order hlij, has the form [19 



2^ 

u 



SO 



TT 



(17) 
(18) 

(19) 
(20) 

(21) 



where — mim2 / {rrii 4-7712) is the reduced mass of the binary, and D is the distance between the source and the 
observer. Q*-' is the quadrupole (or Newtonian) term, P'^-^Q^^ , PQ^^ and P^-^Q''^ are corrections corresponding 
to higher PN orders, PQ^go ^^"^ P^'^Q^so spin-orbit terms [TOl [5]. 

Now we evaluate the polarization states of the detectable gravitational waveform. In order to avoid long 
expressions, we use the components of the vectors N, p, q, v and S in a formal way. After substituting them 
into Eq. and with the use of the transverse-traceless tensor given by Kidder [TH] in the Lense-Thirring 
case we can determine the contributions to the polarization states. To do this we use the decomposition of the 
relative velocity vector ([5]), which suggests a natural structure to /i^^ (see Appendix), and helps us to describe 
the projection in the comoving system. Similarly to /i^j,, the polarization states can be decomposed as 



2m 

Id 



0.5 



(22) 
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where we have used the fact that /U w m in the Lense-Thirring approximation. In this way the expressions for 

the "plus" polarization are: 

^ (Px - qI) + '^v±f{p^Py - qccQy) + vl {pI - ql) 



N 



h 



1.5SO 



= r- - 



(23) 



+ v± 

,.2 



/2Mf 

v'i [rN, + v±Ny]ip'^-q'y) 
1 



M 
2r 



2r ) - 2v±rNy 

,2 „2 



{PxPy - qxqy) 



(24) 



\ 



1 

+ 6 



/2M 

+ Qv±_f h V 

r 



7M^ 



Nl + ^v^r[-^-^ + ?,f']N,Ny 



+ 



(Px - <&) 



{PxPv - qxqy) 



1 

+ 6 



7M 



2vl[-—+ 3f" I N'^ + UvirN^Ny + 6viN^ + vi[-—+3v 



AM 



[(qS)p^ + (pS)9^ 



3r2M v'^M 41M2 



(25) 
(26) 



4r 



12r2 



85Pjy 9^2^ 7_^2 



8r 



8r 



2r2 



3f^ I NiNy 



+ 3..i(?^^-.^)«.«^+.i(^-,.=)«3 

25r27v^ Iv'^M UM'^ 1 



+ r 



Sf^M v'^M 59M2 1 



12r2 



2 -2 



8r 



8r 



3r2 



+ v± 



4r 



/5M 



4r 



2r* 1 Nl + 2t;ir I — - 3r^ ) iV^A^, 



\ 2r 



+ 3t;i — - 2r^ TV^TV^ _ 2t;;[rAf^^ + 



4r 



+ 



4r 



3r2 



2 -2 

V r 



2t)2 Af, 



+ vl 



/29M 



\ 4r 



- 3f^ ) N-^Ny - 3virN,N-^ - vIN^ 



r 
2 



/'3M w2 



(Py - ql) 



[iv^S^{pl - ql) + r[S X (p^p - g^q)]^ - 2t;j_[S x (jj^p - q^ci)]y\ 



For the "cross" polarization we obtain: 



= 2 { f - — 1 p^qx + 2v±r{pxqy + qxPy) + 2v^Pyqy 



(27) 
(28) 

(29) 



/ 4Mf 



+ fx 



r 

/3Af 
V"27 



M 



2r2 



2r A^x - 2v±rNy 



Pxqx 



{Pxqy + qxPy) 



- 2v']_ [tN^, + v±Ny\ pyqy 



(30) 



hi 



( ^ + ^ + 6^ + ^ j + 4«±r- (^-— + 3r2 j AT.AT, 
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ISO 



1.550 



2vi I 3f2 



M 
r 

5M 



Nf. 



IQv^M 29AP 



M 



PxQx 



6v^r [ - — + 2^ ] + 8vi[ -4— + 3f^ 1 N^Ny + l2vlrN^ 



6v±r 



(PxQy + qxPy) 



7M 



2vi + 3P + 12virN^Ny + 6vlN^ + < + 3v 



AM 



Pyly 



[{qS)q., - (pS)p,] 



3PM 2v^M 41M^ 



(31) 
(32) 



2r 



Qj.2 



6rv 



+ v± 



+ v±_ 



f2M 



2P \Nl^v^ 



M 



\ 4r 



4r 

IOPM 2v^M 59M'^ 



6P 



25PM 7v^M 26M^ 



Ar Ar 3P 

PM 7v^M lUP 



v^P N, 



PxQx 



Ar 



Ar 



3ui 



2P N^N'y - 2vlrm + - 



fbM 

f{ 2v^ 

. f5M 

-r 

V 2r 

flAM A /3M 2 



2P ) + 2v^r - 3Pj N^Ny 

A9PM llv^M 32M2 



Ar 



Ar 



3r2 



v'P 



+ 2P ]n;^ + v± 



{Pxqy +qxPy) 

/29M 



\ 2r 



6P N^Ny - 6virN:rNy ~ 2v\N[ 



3 j\r3 

y 



V 2r 



Ny PyQy 



[6v±SzPxqx +r[S x {pxq + qxp)]x - 2v±[S x {pxq + qxp)]y] 



(33) 
(34) 



Eqs. (|23 |) -()34 p are the basics of our main results to determine the time (or an appropriate parameter) 
dependence of the polarization states. After solving the radial equation of motion ((4]) it can be used to solve the 
angular equations ^ and ([7]), and determine the components of the relative velocity vector. Having in hand 
the full description of the motion (up to 1.5 PN order and linear in spin) one can express the components of the 
vectors N, p, q and S in the comoving coordinate system. With the substitution of these vector components into 
Eqs. (|23p - (|34p one can determine the explicit dependence of the polarization states on time (or an appropriate 
parameter). 

To collect all the contributions due to spin effects one has to substitute the decomposition of r, v, r and v±, 
as given in Eqs. (|10p . into the expressions of the polarization states, which results 



[(qS)p, + {pS)q^] 



— [3v±oS:,{pI - ql) + (f)o[S X {p^p - q.j;q)]x - 2u_lo[S x (p^p - feq)]^] 



+ 2 



{r)oir)i + ] {pI ~ qi) + (w±o(?')i + w±i(»')o)(Pa:Py - qxqy) + w±ow±i(Py - 9^) 



S'(35) 



1 [(qS)g, - (pS)p,] 



[Qv^oSzPxqx + ('^)o[S X {pxC\ + qxV)\x - 2t;±o[S X {Px(l + qxP)]y\ 



2{f)o{r)i + Pxqx + {v±Q{r)i + v±i{r)o){pa;qy + qxPy) + 2v±QV±iPyqy 

^0 / 



S , 



(36) 



where (f)i and v_ii denote the coefficients of the terms linear in spin in Eqs. (jlOp . We note that other terms 
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linear in S arise when we consider the evolution of the angles appearing in p and q which are not listed here 
explicitly. 



5 The S = case 

In the S = case the description of the motion can be reformulated. Since the direction of L^r do not change, 
it is useful to choose the z axis of the invariant and the comoving systems to be the same. Formally it means 
that the angle ln is zero, and the separation vector in the invariant system becomes 

cos $ cos vj/ - sin <i> sin \ / cos ($ + 

sin <i> cos * + cos $ sin * =r sin ($ + *) | . (37) 

y V 

It is convenient to introduce the angle T = (f> + "J*. The equation of motion for it is 

t = . (38) 
The radial equation of the motion becomes = (j')q, while the relative velocity vector is 



L 

mr 

/ \ 



rr \ = \ ^ \ ■ (39) 



The X and y axes of the invariant system are again determined by the vector N. Since N — sin 7 in the comoving 
system, the vectors N, p and q become: 

cos T sin 7 \ / sin T \ / — cos T cos 7 

N= I — sin T sin 7 , P= cosT , 1= sin T cos 7 | . (40) 
cos 7 J \ ^ J \ sin7 

Although some of the formal expressions do not change in the spinless limit as compared to the general case, 
one has to substitute 5* = into the different expressions. Also the spin-orbit terms in and hx automatically 
vanish. 



6 The circular orbit case 

A general decomposition of the relative velocity vector is used to describe the circular orbit case: 



7'n + rujin 



(41) 



Using Euler angles the frequency lo has the form 

UJ = cos LN^ + * ■ (42) 

The circular orbit case corresponds to r = and uj = 0. If we use these rules and the method given in jTH] 
for the length of the relative velocity vector we get the following formula 



2 2 
r UJ 



M 
r 



2 

M 



2 (n X m)S 



3/2 



(43) 



Substituting Eq. (^5]) for the components of the relative velocity vector into Eqs. we get the polar- 

ization states of the detectable gravitational waves in the circular orbit limit. Again, we use the decomposition 
of /i+ and /ix as given in Eq. (j22p . The individual terms which correspond to different PN orders and spin 
effects are 



(44) 



7 



and 



1.5SO 



i ISO 
1,1.5 



h 



l.bSO 



1 /M 



3/2 r 



;Ny{pl - ql) + iN^ip^py - q^qy) - Ny{pl ~ q^) 



- ^ — j [{lON^ - 2Nl + 19) {pi - ql) - ?,2N,Ny{p,Py - q^qy) 
{-lANl + ml - 19) [pI - ql)] 



ISO _ 




















+ 


—N. 



M2 



[(qS)p, + (pS)g,] 



5/2 



8 ^ ^ 4 2^ 



101 
'24 



Nyipl 



12 ^ 



15 ..2 

2 y 



Nf, 



19 

T 



Nyipl- ql) 



5/2 



[(2(p2 _ _ (p2 _ ^2)) _ ^^^^g^ 



M 
r 

3/2 



Pylyl 

; + iNxiPxqy + g^Py) - 2NyPyqy] 



1 fM 
3 W 



[{iml - 2iV2 + 19) pxqx - l6NxNy{pxqy + qxPy) 



-UNl+ml~l9) Pyqy] 



M 



-) ^ [(qS)g. - (pS)p.] 



t) 



5/2 



37 



1^,2 101 



-K + l^N^- — ]NyPxq.+ 

y 



12 



65 



15 



149 
12" 



149 



Y^^. + -N'y - — ) N^ip.qy + q.Py) 



12 



- 2Nl + ^ I NyPyqy 



M 



5/2 



= \~] J^[Sz{Pxqx -Pyqy) +PxqzSx] 



(45) 

(46) 
(47) 

NxiPxPy - qxqy) 

(48) 
(49) 



(50) 
(51) 

(52) 
(53) 

) 

(54) 
(55) 



In the circular orbit case one can make the time dependence exphcit since the length of the separation vector 
is constant. The equations of motion can be solved and the polarization states can be given in terms of time 
following the method described in the general case. 

Eqs. (|44p - ([55l) are found to be in agreement with the results given in [12]. To see this one has to reexpress 
Eqs. (4.8) and (4.9) of J^, using the Lense-Thirring approximation and evaluate the projections given in ([1^. 



7 The relation to existing results 

The properties of gravitational waves generated by compact binaries have been studied extensively. Typical 
parameters of the binary, masses, eccentricity and spins, etc. are included in these descriptions. To find the 
relation between our results and the ones known in the literature, we compare the method presented here with 
integrated expressions of the waveforms depending explicitly on time or a chosen parameter. In some cases the 
metric perturbations hl^j, are given formally in terms of the separation vector, the relative velocity and spin 
vectors and N [^. To express the polarization states /i+ and hx one has to introduce suitable coordinate 
systems and integrate the equations of the motion for the angle variables. For this purpose we have used Euler 
angles and the results of P^ . 

Gopakumar and Iyer have obtained the polarization states for eccentric binaries up to 2 PN order. To the 
required accuracy they have included all radiative moments [26j and employed the 2 PN accurate generalized 
quasi-Keplerian parameterization for non-spinning bodies [H]. Here we have used the 1.5 PN expression of 
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h^j,rp and the lowest order spin contributions to the motion of the binary in the Lense-Thirring approximation. 
Higher order corrections to the waveform which arise from an appropriate parameterization of the motion can be 
investigated by substituting the result of the PN expansion of the motion into the quadrupole formula [551 US] ■ 

The effects caused by the rotation of the bodies on the waveform are studied e.g. for circular orbits [20l [27] 
or in the case of binaries with equal masses and one spinning body |28 . 

Most of the results containing explicit time dependent contributions to the waveforms are valid in the 
circular orbit limit. To check our formulae in this case we have used the results of Kidder [19j . where the hl^j, 
tensor is given after the decomposition of the velocity vector, Eq. ([1T|) . As an other check we have integrated 
the equations of the motion in the circular orbit limit. With the use of the method presented above we have 
determined the lowest order contributions to the detectable waveforms. Our expressions reproduce the well 
known result of the quadrupole formalism, namely that the frequency of gravitational waves is twice the orbital 
frequency. It is worth mentioning that the eccentricity of the orbit and the lowest order PN perturbation give 
rise to higher harmonics. 

8 Conclusions and remarks 

We have presented a method to obtain the detectable gravitational wave signals generated by a spinning compact 
binary system in the Lense-Thirring approximation. This is useful for the investigation of the effects caused 
by the rotation of the components in the case of eccentric orbits. We have introduced a comoving coordinate 
system and shown how to express the wave polarizations in this system. Our formal expressions are determined 
in terms of the components of the separation vector, the relative velocity, the spin and the orthonormal triad 
(N, p, q). We have calculated these vector components in terms of the angles describing the motion and the 
relative position of the source and the detector. 

We have clarified all the steps needed to determine the explicit parameter dependence of /i+ and ft-x up 
to 1.5 PN relative order. To do this one can use the generalized true anomaly parameterization of the radial 
motion [23] , which is a natural parameterization of an eccentric orbit in the case of a spinning binary. 

We have investigated the main characteristics of the spinless and the circular orbit limits of our general 
method. We have explored the effects of the spin and the eccentricity of the orbit. The results in the circular 
orbit limit can help one to describe the emitted waves near the innermost stable circular orbit and provide 
initial information for the methods investigating the coalescence of the binary. 

After calculating the time dependence of the detectable gravitational waves following the method given above 
one can construct an experimental receipt to determine the main properties of a compact binary system with the 
detection of the waves it emits. Although recent experiments cannot detect waves emitted by a Lense-Thirring 
system, the method can be useful starting point to investigate the description of the gravitational wave signals 
in more general cases analytically. 
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A The hj^j^ tensor in the Lense-Thirring case 

If we use the Lense-Thirring approximation, and decompose the relative velocity vector according to Eq. ([5|), 
the various contributions of the hl[rp tensor, given in [191, representing different PN orders and effects become 
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